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We present an overview of selected aspects of ultrarelativistic nucleus-nucleus 
collisions, a research program devoted to the study of strongly interacting 
matter at high energy densities and in particular to the characterization of the 
quark-gluon plasma (QGP). The basic features of the phase diagram of nuclear 
matter, as currently understood theoretically, are discussed. The experimental 
program, carried out over a broad energy domain at various accelerators, is 
briefly reviewed, with an emphasis on the global characterization of nucleus- 
nucleus collisions. Two particular aspects are treated in more detail: i) the 
application of statistical models to a phenomenological description of particle 
production and the information it provides on the phase diagram; ii) the 
production of hadrons carrying charm quarks as messengers from the QGP 
phase. 

Go for the messes - that's where the action is. 
S.Weinberg, Nature 426 (2003) 389 



1 Introduction 

Quantum Chromodynamics (QCD), the theory of strong interactions (see [1] 
for a recent review) , predicts a phase transition from a state of hadronic con- 
stituents to a plasma of deconfined quarks and gluons, as the energy density 
exceeds a critical value. The opposite phase transition, from quarks and gluons 
to hadronic matter, took place about 10~^ s after the Big Bang, the primeval 
event which is at the origin of our Universe. The core of the physics program 
of ultrarelativistic nucleus-nucleus collisions research [2] is the study of the 
properties of strongly interacting matter at high energy density, in particular 
its phase diagram and the properties of quark-gluon plasma (QGP) [3]. 

Already in 1951 Pomeranchuk [4] conjectured that a finite hadron size im- 
plies a critical density, ric, above which nuclear matter cannot be in a hadronic 
state. In 1965, Hagedorn [5] inferred that an exponentially growing mass spec- 
trum of hadronic states (observed up to masses of about 1.5 GeV) implies a 
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critical temperature of the order of 200 MeV 2 • 10^^ K). However, the 

elementary building blocks of QCD, the quarks and gluons (carrying an extra 
quantum number called "color") have not been directly observed in experi- 
ments, although their fingerprints have been clearly identified in deep- inelastic 
collisions and jet production. A fundamental property of QCD, the asymptotic 
freedom, unraveled by Gross, Wilczek, and Politzer in 1973 [6], implies that 
the attractive force (coupling) between quarks increases as a function of their 
separation. Moreover, the confinement of quarks (and gluons) inside hadrons 
is another fundamental feature of QCD, although not fully understood yet. 
Cabibbo and Parisi [7] demonstrated already in 1975 that the exponential 
mass spectrum of hadronic states is a feature of any hadronic system which 
undergoes a second order phase transition with critical temperature Tc, since 
thermodynamical quantities exhibit singularities at Tc- This is realized in 
models that include "quark containment" [7], so is in agreement with QCD 
principles. Collins and Perry [8] demonstrated in the same year that asymp- 
totically free QCD is also realized for large densities. It is interesting to note 
that rcf. [7] contains the first sketch of a phase diagram of nuclear matter. 
The term quark-gluon plasma along with initial ideas about the space-time 
picture of hadronic collisions were first introduced by Shuryak [9]. 

2 Theoretical background 

In the recent years, a successful effort to solve the QCD equations numerically 
on a (space-time) lattice has brought deeper insight into the subject of phase 
transition(s) from hadronic to quark-gluon matter [10]. It is not yet clear 
whether the transition is a true singular behaviour of thermodynamic variables 
or just a rapid crossover. 

Fig. 1 shows a map of the phase transition in the coordinates of u,d and 
s quark masses and the chemical potential n [10]. The surface, corresponding 
to a second order phase transition, is the border between the regions of first 
order transition and crossover. While the u and d quarks have small masses 
(of the order of a few MeV), the mass of the s quark is not well known, but 
is likely to be about 150 MeV. In this case (so-called physical values of the 
quark masses, represented by the vertical line in Fig. 1) the transition from 
QGP to a hadron gas is a crossover for small jjb and reaches into the domain 
of the first order for large /i, implying the existence of a critical point [11]. Up 
to now, experimental searches for such a critical point via enhanced event-by- 
event fluctuations have not turned up any evidence [12]. Whether this means 
that all critical fluctuations are effectively damped by the phase transition or 
whether the transition is of first order, is currently an open question. It may 
also imply that the critical point is not reached in the energy range studied 
up to now. 

In any case, the transition to free quarks and gluons is illustrated by the 
sudden increase of the energy density as a function of temperature, shown 
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Fig. 2. Energy density as a function of temperature calculated with lattice QCD 
(taken from ref. [10]). 



in Fig. 2 for two and three degenerate flavors [10]. For the 2-flavor case, the 
transition corresponds to a critical temperature Tc ~170 MeV with critical 
energy density Sc —0.7 GeV, while for the 3-flavor case Tc is smaller by about 
20 MeV. A result for the case of two degenerate flavors and a heavier strange 
quark (physical values) is also included. 

Other features of Fig. 2 can be understood be recalling fundamental results 
of thermodynamics of relativistic gases [13]. The grand partition functions for 
fermions (particles and anti-particles) and bosons are: 
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where gf and gb are the respective degeneracies (degrees of freedom). The 
average energy, particle number and entropy densities and the pressure are: 

_Td{TlnZ) _ 1 (9(rinZ) 9(TlnZ) _ 1 a(rinZ) 

'~v dT "^'v ' ay ' y 

(2) 

Using the thermodynamic relation: e = —P + Ts + jin one can easily estab- 
lish the equation of state (EoS) of an ideal gas: P — e/3. Assuming that the 
hadronic world is composed of pions, gh=^- For three colours and two spin 
values, for quarks and gluons one has gq—12Nf and gg=16, respectively. Nf 
is the number of flavours (the lighter quarks u, d and s are the only ones rele- 
vant). Consequently, at vanishing the chemical potential, the energy densities 
for the hadronic stage and for a gas of free quarks and gluons are, respectively: 

WT' = ^, £,,/T'* = (32-|-21iV/)^. (3) 

For Nf=3, egg/T^=15.6, denoted as the Stefan-Boltzmann limit, Esb, in 
Fig. 2. It is interesting to note that the calculated values are well below the 
values for non-interacting gases, indicating that the QGP is far from an ideal 
gas at temperatures as high as several times Tc- 

An important (and not yet understood) result of lattice QCD calculations 
is that the critical temperatures for deconfinement and for chiral symmetry 
restoration (T^) apparently coincide, although one might expect that > Tc 
[10]. 
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Fig. 3. Schematic phase diagram of nuclear matter (taken from ref. [10]). 
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A simple way to incorporate the two basic properties of QCD, asymp- 
totic freedom and confinement, is achieved in the so-called (MIT) bag model 
[13]. It prohibits quarks ang gluons from existing outside the bag (which can 
be any finite volume) by adding a shift from the physical vacuum into the 
QCD vacuum by an extra term in the partition function of the plasma phase: 
{T\nZ)vac = —BV, where B is the bag constant. It is easy to show that the 
EoS in this case becomes: P = {e — 4B)/3. The phase transition trajectory 
in the T — ^ plane can be constructed by applying the Gibbs criteria for the 
phase transition: 

Ph = Pqg, IJ-h = IJ-qg{= 3M«)) ^ft = ^gg = ^c- (4) 

A sketch of the present understanding [10] of the phase diagram of strongly 
interacting matter is shown in Fig. 3 in the T — n plane, for two light u and 
d quarks and a heavy s quark. The lines mark the borders between the dif- 
ferent phases of hadronic matter. The dots mark the expected position of 
critical points, namely the T — ^ loci beyond which a first order phase transi- 
tion is no longer expected to take place. Ground-state nuclear matter (atomic 
nuclei) corresponds to /io=931 MeV (bomid micleon mass) and r=0 and is 
well modeled as a liquid. The line starting at this point denotes the liquid-gas 
phase boundary which is under study in low energy nucleus-nucleus collisions. 
The region of high temperatures is the part which is being explored in ultra- 
relativistic nucleus-nucleus collisions. The exotic region of low temperatures 
and high densities (high ji) is of relevance to astrophysical phenomena, but is 
rather likely to remain inaccessible to laboratory experiments. 

3 Experimental program and global observables 

By colliding heavy ions at ultrarelativistic energies, one expects to create 
matter under conditions that are sufficient for deconfinement [2]. A series of 
conferences (the so-called "Quark Matter" conferences, see ref. [15] for the 
most recent of those) devoted to the subject has begun in 1980 (Bielefeld). 

The experimental program has started at the CERN's Super Proton Syn- 
chrotron (BPS) and at Brookhaven's Alternating Gradient Synchrotron (AGS) 
in 1985. The AGS program [16], carried out over a period of about 15 years by 
several experiments (E802/864,917E810, E814/877, E864, E895) is essentially 
completed. The SPS program is just being concluded. Compelling evidence 
for the production of a "New State of Matter" , has been produced in cen- 
tral Pb-Pb collisions [17] studied by seven experiments: WA80/98, NA35/49, 
NA38/50/60, NA44 NA45/CERES, WA97/NA57, and NA52. A vigorous re- 
search program, started with the first data taking in 2001, is on-going at the 
Rclativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory 
(BNL) [18] with four experiments, BRAHMS, PHENIX, PHOBOS and STAR. 
The Large Hadron Collider (LHC) will start operating at CERN in 2007 and 
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will provide (in addition to proton beams) heavy ion beams, which will be 

used in the research program of the dedicated ALICE experiment as well as 
by the ATLAS and CMS experiments [19]. A dedicated fixed-target facility 
is planned at Gesellschaft fiir Schwerionenforschung (GSI), expected to be 
operational in 2012 [20]. 

The temporal evolution of a (central) nucleus-nucleus collision at ultra- 
relativistic energies is understood to proceed through the following stages: 
i) liberation of quarks and gluons due to the high energy deposited in the 
overlap region of the two nuclei; ii) equilibration of quarks and gluons; iii) 
crossing of the phase boundary and hadronization; iv) freeze-out. Interest- 
ing experimental information is contained in the study of the distributions of 
(mostly charged) hadrons after freeze-out. Whether any information on the 
phase transition can be gleaned from these investigations will be discussed 
below. Clearly, given the short timescales of a nucleus-nucleus collision and 
the small volume involved (lattice QCD calculations discussed in the previous 
section are for bulk) the reconstruction of the various stages of the collision 
is a difficult task. It is consequently of particular relevance to find experimen- 
tal observables which carry information (preferentially) from one particular 
stage, in particular about the QGP phase. Specific probes of QGP have been 
proposed [21, 22] and are currently being studied experimentally: i) direct 
photons [23]; ii) low- mass dileptons [24]; iii) strangeness [25]; v) charmonium 
suppression [3]; vi) jet-quenching [26]; vii) fluctuations [12, 27]. 

As shown in the next section, the study of hadron multiplicities in a statis- 
tical model is a unique way to provide experimental information on the QCD 
phase diagram [28]. Other global observables, like the distribution of parti- 
cles over momentum space, collective fiow, and the measurements of effective 
source sizes via particle interferometry, have also been studied in detail. In 
particular their energy evolution is of relevance and is briefly examined in the 
following {^snn is the total ccnter-of-mass energy per nuclcon pair). 

In Fig. 4 we present a compilation of experimental data on charged particle 
rapidity density distributions, dNch/dy, and transverse energy rapidity den- 
sity, dET/dy, at midrapidity'^ . The values are for central collisions (average 
value of the number of participant nucleons in the fireball, iVparf=350, which 
roughly corresponds to the 5% most central collisions) in the energy range 
from AGS up to RHIC.^ The continuous lines arc {^ysNN)'^'^ dependences, 
arbitrarily normalized. These power-law dependences describe the measure- 
ments quite well starting from the top AGS energy {^snn —5 GeV). This may 
suggest that some aspects of the imderlying physics are similar over all this 
energy domain. Note that the SPS data (NA49) seem to slightly deviate from 

^ Midrapidity is the rapidity of center-of-mass system; rapidity is defined as y = 
0.51n[(i5 +pz)/{E — Pz)], where pz is the longitudinal component of the particle 
momentum and E is the energy. 

^ A constant Jacobian of 1.1 has been used to convert the dX/dr? data to dX/dy. 
T) = — ln[tan(^/2)] is the pseudo-rapidity {9 is the polar angle of a given particle). 
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Fig. 4. Excitation function 
of global observables in cen- 
tral central nucleus-nucleus col- 
lisions (A'part=350). The exper- 
imental values for particle ra- 
pidity density, dNch/dy [29, 30, 
31, 32, 33, 34, 35] (upper panel) 
and transverse energy rapidity 
density, dEr/dy [30, 36] (lower 
panel) at midrapidity are plot- 
ted as symbols. The lines are 
a power law dependence arbi- 
trarily scaled. The thick hor- 
izontal lines mark the energy 
range of the various accelera- 
tors. The dotted line marks the 
full LHC energy for Pb-Pb col- 
lisions (yijvjv=5.5 TeV). 
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the power-law behaviour. Also, at the lower SPS energies there is an appar- 
ent disagreement between NA49 and NA50/NA60 data. The (^/s]v7v)°'^ de- 
pendences allow for simple, experimentally-based, extrapolations to the LHC 
energy (of course, surprises are eagerly awaited). We note that power-law de- 
pendences are predicted by the (QCD-inspired) saturation model [37], but 
they are steeper (exponent 0.41, in case of dNch/dy). The steep decrease of 
particle multiplicities towards the lower end of the AGS energy range reflects 
mainly the threshold in the overall particle production, but may indicate a 
change in physics as well. The average transverse energy per charged particle 
has a nearly constant value of 0.8-0.9 GeV all the way from top AGS to RHIC 
energies. 

The initial energy density, e, and the net baryon density, ntaryon produced 
in a (central) heavy ion collision can be calculated from the measured trans- 
verse energy (dET/dr/) and net baryon (diVj,_j,/dry) densities, respectively, in 
the so-called " Bjorken-scenario" [21]. This assumes self-similar (Hubble- like) 
homogeneous (hydrodynamical) expansion of the fireball in the longitudinal 
(beam) direction. The resulting relations are: 

^^J_dETdri 1 djV,_fc d?? 

At drj dz ' "-^yon 
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where At is the transverse area of the fireball {At =154 fm^ for a head-on 
Au-Au collision). In the above equations the only unknown parameter is the 
formation time (the time for establishing the equilibrium), r (djj/dz = 1/t), 
which is usually taken to be 1 fm/c, although it is expected to decrease as a 
function of the energy. In this sense, the values obtained using Eq. 5 arc con- 
servative estimates for most of the energy range spanned by the experiments. 



Table 1. Measured and deduced quantities at AGS, SPS, and RHIC for central 



Machine 


AGS SPS RHIC LHC 


^snn (<-ieVj 


4.9 17.3 200 5500 


dEr/dri (GeV) 
dN,_-Jdv 


192 363 625 1800 ? 
170 100 25 -0 ? 


£ (GeV/fm^) 

IT'baryon (fni ) 


1.2 2.4 4.1 11.6 ? 
1.1 0.65 0.17 ? 



The maximum nucleon-nucleon center-of-mass energy {^/snn) and the 
corresponding measured and calculated (using Eq. 5) energy and baryon den- 
sities arc listed in Table 1 for the various accelerator regimes. The results for 
LHC are extrapolations based on the (^/snn)^-^ dependence discussed above. 

The achieved densities are obviously very much larger than those inside 
a normal Pb nucleus: £o = 0.15 GeV/fm^ and no = 0.16/fm^. The estimates 
of the energy densities are for all the energy range above the critical energy 
density a.t /j-b = (ec — 0.7 GeV/fm^), indicating that the conditions for the 
QGP formation likely have been achieved in the experiments. 

In Fig. 5 we present the excitation function of elliptic flow [38] for semi- 
central collisions (for which this observable has a maximum as a function 
of centrality). This observable is characterized by the second order Fourier 
coefficient V2 = (cos((/))), where (j) is the azimuthal angle with respect to 
the reaction plane (defined by the impact parameter vector and the beam 
direction). It reflects the initial geometry of the overlap region and its pressure 
gradients. A transition from out-of-plane (also called "squeeze-out", V2 <0) to 
in-plane (v2 >0) preferential particle emission is seen in the energy domain of 
the AGS. This is the result of a combined effect of the violence of the expansion 
and of the shadowing of the spectator matter, which, at these energies is still 
present in the vicinity of the fireball. A striking correlation of this transition 
with the sharp increase of particle multiplicity at midrapidity (seen in Fig. 4) 
is evident. From top AGS energy up to RHIC the V2 values increase steadily 
and are also well described by a log(ys]v7v) dependence. An early observation 
at RHIC was that the V2 values are reaching the hydrodynamical [39] limits 
which is an indication of an early equilibration of the fireball. It is likely 
that QGP is the only way to achieve such a fast equilibration. It is thus an 
interesting question whether elliptic flow at LHC will follow the log(^sjvjv) 
trend or will flatten at the RHIC values. 
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Fig. 5. Excitation function of elliptic flow. Protons are considered up to SPS energies 
and all charged particles at RHIC. The line is a log(-^/sjvjv) dependence arbitrarily 
normalized. 
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Fig. 6. Excitation function of the freeze-out volume extracted from pion HBT 
correlations. The line is a log(y^sjvjv) dependence arbitrarily normalized. 



In Fig. 6 we show the energy dependence of the volume of the fireball as 
extracted from pion Hanbury Brown-Twiss (HBT) correlations [40]. Again, 
a strikingly different behavior is seen at the lowest energies compared to top 
AGS and above, for which a log(y^s]v7v) dependence describe the measure- 
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ments [41] well. This non-monotonic behavior can be understood [42] quan- 
titatively as the result of an universal pion freeze-out at a critical mean free 
path Xf —1 fm, independent of energy. It is worth mentioning that the smooth 
evolution of the source size in the energy range of top AGS to RHIC is an 
indication that no first order phase transition, associated with supercooling 
and explosive expansion, is visible in hadronic observables in this energy do- 
main. Also, the measured source sizes at RHIC are well below hydrodynamic 
predictions. 

4 Particle yields and their statistical description 

The equilibrium behavior of thermodynamical observables can be evaluated 

as an average over statistical ensembles. The equilibrium distribution is thus 
obtained by an average over all accessible phase space. Furthermore, the en- 
semble corresponding to thermodynamic equilibrium is that for which the 
phase space density is uniform over the accessible phase space. In this sense, 
filling the accessible phase space uniformly is both a necessary and a sufficient 
condition for equilibrium. We restrict ourselves here on the basic features and 
essential results of the statistical model approach. A complete survey of the 
assumptions and results, as well as of the relevant references, is available in 



The basic quantity required to compute the thermal composition of par- 
ticle yields measured in heavy ion collisions is the partition function Z{T, V). 
In the grand canonical (GC) ensemble, for particle i of strangeness Si, 
baryon number Bi, electric charge Qi and spin-isospin degeneracy factor 
Qi = (2Jj + l)(2/j -I- 1), the partition function is: 



with for fermions (like baryons, made of 3 quarks) and (-) for bosons (like 
mesons, made of quark-antiquark pairs). Note that the partition functions in- 
troduced in Eq. 1 are for massless particles, for which the analytic integration 
of Eq. 6 can be performed. The particle density is: 



T is the temperature and Ei = y/p'^ + mf is the total energy, jii = + 
lis Si -\- fii-ihi is the chemical potential, with ^ib, ^^^d ^iq the chemi- 
cal potentials related to baryon number, strangeness and electric charge;, rc;- 
spectively, which ensure the conservation (on average) the respective quan- 
tum numbers: i) baryon number: V ^ - UiBi = Z + N; ii) strangeness: 
V^-UiSi = 0; iii) charge: V^^riil^i = . This leaves T and the bary- 
ochemical potential fij, as the only parameters of the model. In practice, how- 



ref. [28]. 
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ever, the volume determination may be subject to uncertainties due to incom- 
plete stopping of the colliding nuc;lci. Due to this reason, the most convenient 
way to compare with measurements is to use particle ratios. 

The interaction of hadrons and resonances is usually included by imple- 
menting a hard core repulsion of Van dcr Waals typo via an excluded volume 
correction. This is implemented in an iterative procedure according to: 

p«^^'-(r,M) = p*'^-9'^^(T,A); ^^ = ^Ji-Ve,genP'''''■{T,^l) (8) 

where Veigen is calculated for a radius of 0.3 fm, considered identical for all 
particles. 

The grand canonical ensemble is of course the simplest realization of a 
statistical approach and is suited for large systems, with large number of 
produced particles. However, for small systems (or peripheral nucleus-nucleus 
collisions) and for low energies in case of strangeness production, a canoni- 
cal ensemble (C) treatment is mandatory. It leads to severe phase space re- 
duction for particle production (so-called "canonical suppression"). Within 
this approach, particle production in e+e" collisions has been successfully de- 
scribed, albeit with an additional heuristic strangeness suppression factor. It 
has been shown that the density of particle i with strangeness S calculated in 
the canonical approach, np, is related to the grand canonical value, np*^, as; 
nf = nf^Fs, with Fs = Is{x)/ Iq{x). The argument of the Bessel function of 
order S is the total yield of strange and antistrange particles. For central Pb- 
Pb (Au-Au) collisions, the canonical suppression is negligible for all strange 
particle species already for the highest AGS energy {yfsNN —5 GeV) but is 
sizeable for the lowest energy considered in the following, ^Jsnn =2.7 GeV 
(corresponding to the beam energy of 2 GeV/n), for which Fi ~2, F2 —8. 

In Fig. 7 we present the result of a thermal fit of the measured parti- 
cle ratios for Pb-Pb collisions at 158 GeV/nucleon beam energy. The values 
r=170±5 MeV and ^f,=255±10 MeV arc the free parameters. The reduced 

(excluding (j) and d) is 2.0, of which the largest contribution comes from 
the ratios vI/tt, A/h~ and yl/K°, possibly due to weak decays feeding. 

The thermal fits of particle ratios for the RHIC energies (ys]vjv=130 
and 200 GeV) are shown in Fig. 8. The obtained values for (T,/Xb) are 
(174±7,46±5) MeV and (177±7,29±6) MeV, respectively with reduced 
values of 0.8 and 1.1. 

We mention here that the measured enhancement of strange hyperons {A, 
E, n) at SPS in central Pb-Pb collisions with respect to pBe and pPb (a factor 
of 20 enhancement in case of fi) can be understood quantitatively not as an 
enhancement in central Pb-Pb but as a suppression in pBe/pPb with respect 
to central Pb-Pb. It is also important to note that, at RHIC, the transverse 
momcntiim spectra can be well described in a thermal approach, with two 
additional (size) parameters [44]. At AGS, the measured yields of light nuclei 
{A <7) are well explained by the thermal model [43]. 

T and /it were determined for other energies (SPS at 40 GeV/n, AGS at 
10.8 GeV/n and for 1 GeV/n Au-Au collisions at SIS) with a similar fitting 
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Fig. 7. Fit of particle ratios for Pb-Pb collisions at SPS (158 GeV/c). The mea- 
surements are the symbols, the thermal fit values are the lines. 
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Fig. 8. Fit of particle ratios for Au-Au collisions at RHIC. The measurements are 
the symbols, the thermal model values are the lines. 



procedure, although using in most cases fewer available measured ratios [28]. 
The resulting values are shown in a phase diagram of hadronic matter [43] 
in Fig. 9, together with calculations of freeze-out trajectories or a hadron gas 
at constant energy density and at constant baryon density. This latter case, 
corresponding to nb=0.12 fm"'^, does reproduce well the freeze-out points 
extracted from the data. Another observation is that the freeze-out points 
lie on a curve corresponding to an average energy (E) per average number 
of hadrons (N) of approximately 1 GeV. We have mentioned above that an 
universal pion freeze-out corresponding to a mean free path of about 1 fm has 
been derived from HBT source size measurements [42]. 
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Fig. 9. The phase diagram of 
nuclear matter in the T — /ij, 
plane. The dots represent the 
extracted values from thermal 
fits to measured particle ratios. 
The trajectories of freeze-out 
for a hadron gas at constant en- 
ergy density (e=500 MeV/fm^) 
and at constant baryon den- 
sity (ni,=0.12 fm~^) are shown 
by the dotted and dash-dotted 
lines, respectively. The phase 
boundary from lattice QCD 
(LQCD) calculations is shown 
with continuous lines. The 
open squares indicate the criti- 
cal point with two difi'erent in- 
puts for calculations [11]. 
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An important observation about the phase diagram is that, for the top 
SPS energy and above, the thermal parameters are (implying hadron yields 
frozen) at the phase bomidary, as known from lattice QCD calculations [11]. A 
natural question though is how is equilibrium achieved? Considerations about 
collisional rates and timescales of the hadronic fireball expansion [45] imply 
that at SPS and RHIC the equilibrium cannot be established in the hadronic 
medium and that it is the phase transition which drives the particles densities 
and ensures chemical equilibrium. 

In a recent paper [46] many body collisions near Tc were investigated as a 
possible mechanism for the equilibration. There it is argued that because of 
the rapid density change near a phase transition such multi-particle collisions 
provide a natural explanation for the observation of chemical equilibration at 
RHIC energies and lead to T = Tc. to within an accuracy of a few MeV. Any 
scenario with T substantially smaller than Tc would require that either multi- 
particle interactions dominate even much below Tc or that the two-particle 
cross sections are larger than in the vacuum by a high factor. Both of the 
latter hypothesis seem unlikely in view of the rapid density decrease. The 
critical temperature determined from RHIC for T ^ Tc coincides well with 
lattice estimates [10] for /i = 0, as discussed above. The same arguments as 
discussed here for RHIC energy also hold for SPS energies: it is likely that also 
there the phase transition drives the particle densities and ensures chemical 
equilibration. 

We note that thermal models have also been used [47] to describe hadron 
production in e"*"e^ and hadron-hadron collisions, leading to temperature pa- 
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rameters close to 170 MeV. This suggests that hadronization itself can be 
seen as a prethermalization process. However, to account for the strangeness 
undersaturation in such collisions, multi-strange baryons can only be repro- 
duced by introducing a strangeness suppression factor of about 0.5, leading 
to a factor of 8 suppression of fl baryons. This non-equilibrium feature, also 
visible in the momentum distributions of the produced particles, is most likely 
due to the "absence" of multi-particle scattering since the system is not in a 
high density phase due to a phase transition. 
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Fig. 10. Energy dependence of the thermal parameters T and The symbols 
are the values extracted from experimental data, the lines are parametrizations (see 
text). 

The energy dependence of the extracted T and fij, values is presented 
in Fig. 10. The lines are parametrizations that allow for extrapolating the 
parameters up to the LHC energy. For /Zb the following parametrization has 
been used [28]: 

A^b = 1270[MeV]/(l + V^[GeV]/4.3), (9) 

while T has been described with a Fermi-like function. For both cases the 
parametrizations describe well the extracted values over all the energy range. 

In Fig. 1 1 we present excitation functions for a selection of strange particle 
yields over the whole energy range from lowest AGS to LHC energy. The ex- 
perimental ratios /n^ , A/n^ and S~ /A are calculated from measurements 
of absolute yields of 7r± [29, 48, 31, 32, 49, 50], [51, 31, 32, 49, 50], A [52] 
and S~ [53]. The errors reflect the systematic uncertainties. These ratios are 
compared to thermal model calculations employing the parametrizations of T 
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Fig. 11. Excitation function for strange particle production, _ff^ and A yields rela- 
tive to pions, S~ relative to A. All the measured data (symbols) are for midrapidity, 
with the exception of A and S~ at AGS, for which only 47r yields are available. The 
lines are thermal model calculations for three cases of weak decay reconstruction 
efficiencies (see text). 

and fib of Fig. 10. In case of the calculations the contribution (mainly impor- 
tant for pion yields) of down-feeding from resonances (via their weak decays) 
is taken into account in three different cases, assuming that none, 50% or all 
of the weak decays contribute to the yields. As one can see, the effect is sig- 
nificant, implying that it is very important that the experimental conditions 
(vertex cuts for selecting particles) for extracting the yields are well specified 
and taken into account in the model calculations. In a way, the extremes in 
weak decays reconstruction fraction shows the range of systematic uncertain- 
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ties that can arise in the comparison of model results with experimental data, 
if experimental information on feeding is ignored (or not known). 

Given the accuracy of the description of multi-particle ratios presented 
above, it is not surprising that overall the model does reproduce the experi- 
mental values rather well up to RHIC energies. The observed discrepancies can 
be explained by the constant temperature (170 MeV) used for these calcula- 
tions, which is not identical (although close) with the temperatures extracted 
from fits of multiparticle ratios shown above. An apparent disagreement be- 
tween measurements and the model calculations is seen concerning the energy 
dependence of the j-K^ and yl/7r+ ratios at SPS energies. The origin of the 
rather narrow structure in the data is currently much debated [54]. We note 
that transport models also cannot reproduce the i^+/7r+ ratio [55]. 

The four ratios presented in Fig. 11 have a very different dependence on 
energy, which reflects the evolution of the fireball at freeze-out, dominated 
by the initial nucleons at low energies and by the newly created particles 
at RHIC and beyond. At LHC, it is expected that the fireball will consist 
exclusively of created particles. The steep variation of the ratios at the lowest 
energies reflects the close threshold for strangeness production. The canonical 
suppression plays an important role as well. 




2 3 
10 10 10 

Vs,N (GeV) 

With the T and \ih values fixed by the fits to the measured particle ratios 
over a broad energy range, the thermal model has a good predictive power 
for all possible particles that can be formed at freeze-out. As an example, 
predictions for thermal yields of K~ clusters [56] relative to A hyperon are 
shown in Fig. 12 in comparison with the ratio j A. Such exotic K~ bound 
states have been predicted to form due to the strongly attractive K~ poten- 
tial within nuclear matter [56] , but are not yet observed experimentally. The 
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yield of single- -ftT" systems have large values, significantly above S~ yields, at 
low energies and exhibit a pronounced decrease as a function of energy. The 
energy dependence of double-X" systems exhibits a broad maximum around 
^snn —6 GeV, a region which will be covered by the future GSI accelerator 
[20]. 

In closing this section, we note that an open question remains concerning 
statistical model description of strongly decaying resonances (like p meson 
and A baryon). Their yields are strongly underestimated by the calculations 
[58, 28]. 



5 Charmonium and charmed hadrons 

The importance of the so-called hard probes, among which the creation of 
heavy-quarks (c and b) have a prominent place, stems from the fact that they 
are exclusively created in primary hard collisions. Consequently, they are ideal 
messengers of the early stage (QGP phase) of the collision. In particular the 
J/ijj meson, which is a bound state of c and c quarks, was predicted to melt in 
the quark-giuon plasma [3] , thus providing a clear signature of its existence. 
Although recent theoretical investigations based on lattice QCD cast doubt 
on the melting at T < l.bTc [57], there is continued interest in quarkonia as 
probes of the QGP. 
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Fig. 13. Energy dependence 
of the total charm cross sec- 
tion in elementary (pp) colli- 
sions. The measurements per- 
formed with nucleus-nucleus 
experiments are compared to 
NLO pQCD calculations [62] 
for the integral and rapidity 
density cross section. 



The production mechanisms of open charm (D mesons) and open beauty 
(B mesons) in elementary collisions can be well described by perturbative 
QCD (pQCD) calculations. For instance, data on charmed meson production 
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over a broad energy range was found [59] to be in good agreement with cal- 
culations using the PYTHIA code (in leading order approximation, so that a 
scale factor of 5 has been used in ref. [59] to approximate the next-to-leading 
order, NLO). Available experimental data on quarkonia production in pp col- 
lisions {J /ip data for -^s below 100 GcV and the T family data up to Tcvatron 
energy, -^5=1. 8 TeV) have been successfully compared to pQCD calculations 
[60]. Recent measurements of J/V' in pp collisions at RHIC are well described 
by (tuned) pQCD calculations, together with the measurements available at 
lower energies [61]. 

NLO pQCD calculations for total charm cross section in elementary colli- 
sions show clearly that the results depend significantly on the choice of several 
parameters, like the parton distribution function (PDF), charm quark mass 
(mc) and renormalization and factorization constants, hr and jip [62]. This 
dependence is the bigger the larger the energy, so it is most crucial for LHC 
energies. A comparison of these calculations with data is presented in Fig. 13. 
With this choice of parameters, the calculations somewhat underpredict the 
measured values. Note that all the measurements are indirect: at SPS the cross 
section was estimated from the measured Drell-Yan cross section in pp [63], 
while at RHIC it was extracted from the charm contribution to the single- 
electron spectra measured in Au-Au collisions [64]. 

In nucleus- nucleus collisions, the J/^ production at SPS is well measured 
by the NA50 collaboration [65] (see also ref. [66] for an in depth discussion). 
The measured V'VV' ratio [67] is independent of energy and is in p-A collisions 
the same as in pp. This ratio is decreasing as a function of centrality in Pb- 
Pb collisions, as seen in Fig. 14 [67, 68], and reaches, for central collisions, 
a value expected for a thermal ensemble at T ~170 MeV [69, 68]. For an 
interpretation of this result, see below. 
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Fig. 14. Centrality dependence of the ratio / {J / tl)) (including branching ratios 
into At+M") at SPS. 
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At RHIC, the recent measurements of J/ip [70] are hampered by very poor 
statistics, but high quahty data arc expected in the near future. First rcsuhs on 
open charm production at RHIC in d-Au colhsions have just been announced 
[71]. In the near future, open charm cross sections will be extracted from 
recently-completed measurements in In-In collisions by NA60 [34]. At LHC, 
there are good prospects to measure a complete set of charm and bottom 
particles [72], in particular with the ALICE experiment. 

Below we discuss the QGP fingerprints as could be unraveled through 
the model of statistical hadronization of charm quarks [68]. A kinetic model 
description of J/ij} production has been independently developed [73]. It is 
equivalent from the point of view of the physical assumptions with the model 
discussed here, but differs in its numerical realization. Other approaches to 
statistical hadronization exist [74, 75], which differ from the model discussed 
here, but mostly in terms of inputs, while the outcome is qualitatively similar. 
The statistical hadronization model (we follow here the outline of ref. [76]) 
assumes that all charm quarks are produced in primary hard collisions and 
equilibrate^ in the quark-gluon plasma. An important corollary of this as- 
sumption is that no J/ip mesons are preformed in the QGP, implying that 
the dissociation of J/V' in QGP [3] is complete. As noted above, recent lat- 
tice QCD calculations show that the J/t/" mesons may not be dissociated in a 
deconfined medium below about 1.5Tc [57]. However, it is possible that, from 
SPS energy on, the initial temperature achieved in the collision exceeds this 
value. 

The question of charm equilibration is a difficult one, but needs to be 
addressed. The cross sections for production of charmed hadrons are much 
too small [77] to allow for their chemical equilibration in a hadronic gas. But 
how can the apparently "thermal" values of the ratio ip' /iphe reconciled with 
this finding ? We assume that all charm quarks are produced in initial hard 
collisions, but that open and hidden charm hadrons are formed at chemical 
freeze-out according to statistical laws. Consistent with the fact that, at the 
top SPS energy and beyond, chemical freeze-out appears to be at the phase 
boundary (sec previous section), the model implies that a QGP phase was 
a stage in the evolution of the fireball. The analysis of J ftp spectra at SPS 
[78] lends further support to the statistical hadronization picture where J/ip 
decouples at chemical freeze-out. A recent analysis of single-electron spectra 
at RHIC [79] also strenghtens the case for an early thermalization of heavy 
quarks. However, in that analysis it was pointed out that both the hydrody- 
namical approach and PYTHIA reproduce the measured single-electron spec- 
tra, although the two approaches are different in detail at low pt and differ 
manifestly at high pt [pt mass of charm quark). Another theoretical anal- 
ysis [80] indicates though that charm quarks might not thcrmalize quickly 
because of their large mass. All of this emphasizes the need to have high- 

This implies thermal, but not chemical equilibrium for charm quarks. 
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precision direct measurements of open charm, which could impose constraints 
on different interpretations. 

In statistical models charm production needs to be treated within the 
framework of canonical thermodynamics [28] . Thus, the charm balance equa- 
tion required during hadronization is expressed as: 

Here iV^g'' is the number of directly produced cc pairs and /„ are modified 
Bessel functions. In the fireball of volume V the total number of open A"*^ = 
n*^F and hidden N*^ — n'-J^V charm hadrons are computed from their grand- 
canonical densities n*^ and respectively. The densities of different particle 
species in the grand canonical ensemble are calculated following the statistical 
model [28] introduced in the previous section. All known charmed mesons and 
hyperons and their decays are included in the calculations. 

The balance equation (10) defines a fugacity parameter Qc that accounts for 
deviations of charm multiplicity from the value that is expected in complete 
chemical equilibrium. The yield of open charm mesons and hyperons i and of 
charmonia j is obtained from: 



The above model for charm production and hadronization can be only 
used if the number of participating nucleons Npart is sufficiently large. Taking 
into account the measured dependence of the relative yield of ijj' to J/ij} on 
centrality in Pb-Pb collisions at SPS energy, seen in Fig. 14, the model appears 
appropriate for Npart >100, for which the ratio approaches the thermal value 
[69, 68]. 

To calculate the yields of open and hidden charm hadrons for a given 
centrality and collision energy one needs to fix a set of parameters in Eq.(lO) 
and (11): 

i) A constant temperature of 170 MeV and a baryonic chemical poten- 
tial fjLh according to the parametrization (9) are used for our calculations 
(see Fig 10). These thermal parameters are consistent with those required 
to describe experimental data on different hadron yields for SPS and RHIC 
energies. 

ii) The volume of the fireball. We focus on rapidity density calcula- 
tions which are of relevance for the colliders, so in this case the volume 
corresponds to a slice of one rapidity unit at midrapidity, VAy=i- It is ob- 
tained from the charged particle rapidity density ANchl'^Hi via the relation 
dNch/dy = n*^V/ij^=i, where n*^ is the charged particle density computed 
within the thermal model. The charged particle rapidity densities (and total 
yields in case; of SPS, for which we calculate Att yields for a direct compari- 
son to experimental data) are taken from experiments at SPS and RHIC and 
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extrapolated to LHC energy (as seen in Fig. 4). Central collisions correspond 
to 7Vpart=350. For the centrality dependences we assume that the volume of 
the fireball is proportional to Npart- 



Table 2. Input {dNch/dy and AN^"^ /dy) and output {VAy=i and Qc) parameters for 
model calculations at top SPS, RHIC and LHC for central collisions {Npart=i^Q)- 



^SNN (GeV) 


17.3 200 5500 


dA^,„/(f'/ 

dN'J^' /dtj 


4.30 730 2000 
().()() 1 i.!)2 16.8 


VAy=i (hn^) 

gc 


861 1663 4564 
1.86 8.33 23.2 



iii) The yield of open charm dN^^'^/dy at midrapidity (or in full volume) 

is taken from NLO pQCD calculations for pp collisions [62] and scaled to 
nucleus-nucleus collision via the nuclear overlap function, Taa [81]. For a 
given centrality: 

diV^ , _ dajpp^cc) ^ 

y^part) = J-AA(,J>lpart)- (1^) 

The pQCD calculations with the MRST HO PDF are used here. 

The input values dNch/dy and dN^^^ /dy and the corresponding volume 
at midrapidity and enhancement factor are summarized in Table 2 for model 
calculations for different collision energies. 

We first compare predictions of the model to 47r-integrated J/iJj data at 
the SPS measured by NA50 collaboration [65, 82]. For the fireball total vol- 
ume 1^=3070 fm^ (for Nch = 1533) the total yield of thermal open charm 
pairs is iV*^=0.98. This is to be contrasted with N^^=0.137 from NLO calcu- 
lations [62], leading to a value of gc=0-78. Although gc is here close to unity, 
this obviously does not indicate that charm production appears at chemical 
equilibrium, as the suppression factor is a strongly varying function of the col- 
lision energy. We have already indicated that, within the time scales available 
in heavy ion collisions, the chemical equilibration of charm is very unlikely 
both in confined and deconfined media. 

In Fig. 15 we show the comparison between the results of our model and 
NA50 data for two different values of N^^": from NLO calculations [62] and 
scaled up by a factor of 2.8. Using the NLO cross sections for charm production 
scaled by the nuclear overlap function, the model understimates the measured 
yield. To explain the overall magnitude of the data, we need to increase the 
iV^g"" yield by a factor of 2.8 as compared to NLO calculations. We mention in 
this context that the observed [63] enhancement of the di-muon yield at inter- 
mediate masses has been interpreted as a possible indication for an anomalous 
increase of the charm production cross section. A third calculation (resulting 
in the dash-dotted line in Fig. 15) is using the NLO cross section scaled-up 
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Fig. 15. The centrality depen- 
dence of J/ip production at SPS. 
Model predictions are compared 
to 47r-integrated NA50 data [65, 
82]. Two curves for tiie model 
correspond to the values of N^^^ 
from NLO calculations (dashed 
line) and scaled up by a factor of 
2.8 (continuous line). The dash- 
dotted curve is obtained when 
considering the possible NA50 
A'part-dependent charm enhance- 
ment over their extracted pp cross 
section [63] (see text). 
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by 1.6, which is the ratio of the open charm cross section estimated by NA50 
for pp collisions at 450 GeV/c [63] and the present NLO values. For this case 
the Npart scaling is not the overlap function, but is taken according to the 
measured di-muon enhancement as a function of Npart [63] . The resulting J/ip 
yields from the statistical model are on average in agreement with the data, 
albeit with a flatter centrality dependence than by using the nuclear over- 
lap function. Thus our charm enhancement factor of 2.8 needed to explain 
the J/ip data is very similar to the factor needed to explain the intermediate 
mass dilepton enhancement assuming that it arises exclusively from charm 
enhancement [63]. We note, however, that other plausible explanations exist 
of the observed enhancement in terms of thermal radiation [83] . 

We turn now to discuss our model predictions for charmonia and open 
charm production at collider energies and compare them with the results ob- 
tained at SPS. Notice that from now on we focus on rapidity densities, which 
are the relevant observables at the colliders. In Table 3 we summarize the 
yields for a selection of hadrons with open and hidden charm. All predicted 
yields increase strongly with beam energy, reflecting the increasing charm 
cross section and the concomitant importance of statistical recombination. 
Also, ratios of open charm hadrons evolve with increasing energy, reflect- 
ing the corresponding decrease in the charm chemical potential. Very recent 
measurements of open charm in d-Au collisions at RHIC [71] yield the ra- 
tio {D*++B*-)/{J:)°+I)°) of 0.40±0.09, which is in a good agreement to the 
model prediction of 0.42. 

Model predictions for the centrality dependence of J/ip and D+ rapid- 
ity densities normalized to Npart are shown in Fig. 16. The results for J/ip 
mesons exhibit, in addition to the dramatic change in magnitude, a striking 
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Table 3. Mid-rapidity densities for open and hidden charm hadrons, calculated for 
central collisions (iVpart=350) at SPS, RHIC and LHC. 
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Fig. 16. Centrality dependence 
of rapidity densities of J/ip (up- 
per panel) and (lower panel) 
mesons per Npart at SPS, RHIC 
and LHC. Note the scale factors 
for RHIC and SPS energies. The 
dashed lines in the lower panel 
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represent A'^^^^ dependences nor- 
malized for Npart~350. 
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change in the shape of the centrahty dependence. In terms of the model this 
change is a consequence of the transition from the canonical to the grand- 
canonical regime. For D"'"-mesons, the expected approximate scaling of the 
ratio D~^/ Npart oc A^pa^t (dashed lines in Fig. 16) is only roughly fulfilled due 

to departures of the nuclear overlap function from the simple N^j^f. depen- 
dence. 
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The results summarized in Table 3 and shown in Fig. 16 obviously depend 
on two input parameters, dNch/dy and dN^^'^/dy. For LHC energy, neither 
one of these parameters is well known. An increase of charged particle multi- 
plicities by up to a factor of three beyond our "nominal" value dA^c/i/dy=2000 
for central collisions is conceivable. However, due to quite large uncertainties 
on the amount of shadowing at LHC energy, these results may be still modi- 
fied. The yield of dN^^^ /dy is also not well known at LHC energy. Although 
these uncertainties affect considerably the magnitude of the predicted yields, 
their centrality dependence remains qualitatively unchanged: the yields per 
participant are increasing functions of Npart- We also note here that, while de- 
tailed predictions differ significantly, qualitatively similar results (see ref. [72] ) 
have been obtained for a kinetic model study of J /ip production at the LHC. 
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Fig. 17. Centrality dependence 
of rapidity densities of J/ip 
mesons at RHIC. Upper panel: 
sensitivity to dNch/dy; middle 
panel: sensitivity to AN^^^ /Ay; 
lower panel: sensitivity to T. The 
calculations are represented by 
lines. The dots are experimental 
data from the PHENIX collabo- 
ration [70]. Note that the point 
for the central collisions is the up- 
per limit extracted by PHENIX 
for 90% C.L. [70]. 



In Fig. 17 we present the predicted centrality dependence of the J /i}) rapid- 
ity density normalized to Npart for RHIC energy (y^s]vjv=200 GeV). The three 
panels show its sensitivity on dNch/dy, dN^^^ /dy, and (freeze-out) tempera- 
ture T. The calculations are compared to experimental results of the PHENIX 
Collaboration [70] . The experimental data have been rescaled according to our 
procedure to calculate Npart and the number of binary collisions, Ncou- Within 
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the still large experimental error bars, the measurements agree with our model 
predictions. In Fig. 17 only the statistical errors of the mid-ccntral data point 
are plotted. The systematic errors are also large [70]. A stringent test of the 
present model can only be made when high statistics J/ij} data are available. 

Wc turn now to a more detailed discussion of the sensitivity of our cal- 
culations to the various input parameters as quantified in Fig. 17. First we 
consider the influence of a 10% variation of dNch/dy on the centrality depen- 
dence of J/tp yield. Note that the total experimental uncertainty of dNch/drj 
(which is for the moment the measured observable for most experiments) at 
RHIC is below 10%. The sensitivity on the dNch/dy values stems from the 
volume into which the (fixed) initial number of charm quarks is distributed. 
The smaller the particle multiplicities and thus also the fireball volume, the 
more probable it is for charm quarks and antiquarks to combine and form 
quarkonia. That is why one sees, in the top panel in Fig. 17, that the J/ip 
yield is increasing with decreasing charge particle multiplicity. 

The sensitivity of the predicted J/tp yields on dA^*'"/dy is also straight- 
forward. The larger this number is in a fixed volume the larger is the yield 
of charmed hadrons. In case of charmonia the dependence on dN^i^ /dy is 
non-linear due to their double charm quark content, as reflected by the factor 
Qc in equation (11). To illustrate the sensitivity of the model predictions on 
dN^g^ /dy, we exhibit the results of a 20% variation with respect to the value 
given in Table 2. The open charm cross section is not yet measured at RHIC. 
However, some indirect measurements can be well reproduced, within the ex- 
perimental errors, by PYTHIA calculations using a p-p charm cross sections 
scaled with the number of collisions Ncoii of 650 /xb [64]. The corresponding 
value at ^/sn^=130 GeV is 330 iih [64]. For comparison, the NLO pQCD 
values we are using are 390 and 235 /zb, respectively. Despite the still large 
experimental uncertainties, this discrepancy needs to be understood. We note 
that, dependent on the input parameters used in the NLO calculations [62], 
possible variations of the open charm production cross section for the RHIC 
energy are of the order of ±20%. In terms of our model this variation cor- 
responds to about a ±30%) change of the J/ijj yield, which is also centrality 
dependent (see middle panel in Fig. 17). If we use the PHENIX p-p cross 
section of 650 /ub, the calculated yield is a factor 2.5 larger for Npart=350 
and increases somewhat stronger with centrality. As apparent in Fig. 17, the 
predictive power of this model, or of any similar model, relies heavily on the 
accurate knowledge of the charm production cross section. A simultaneous 
description of the centrality dependence of open charm together with J/ip 
production is, in this respect, mandatory to test the concept of the statistical 
origin of open and hidden charm hadrons in heavy ion collisions at relativistic 
energies. 

The apparent weak dependence of J/tp yield on freeze-out temperature, 
seen in Fig. 17, may be surprising. In our model this result is a consequence of 

the charm balance equation (11). The temperature variation leads, obviously, 
to a different number of thermally produced charmed hadrons, but this is 
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compensated by the Qc factor. The approximate temperature dependence of 
Qc and the J/ip yield are: 

g,(T)~l/<~e^, A^,//^(r)=g2^,*)^~e""°"-'"''* (13) 

As a result of the small mass difference in the exponent the J/ tp yield exhibits 
only a weak sensitivity on T. This is in contrast to the purely thermal case 
where the yield scales with exp(— m j/^^/T). The only exception is the ratio 
tp'/J/fpj which is obviously identical in the statistical hadronization scenario 
and in the thermal model and coincides, for T~170 MeV, with the measured 
value at SPS (see Fig. 14). 
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Fig. 18. Centrality dependence 
of rapidity densities of J/xp 
mesons at RHIC for different ra- 
pidity window sizes. The lines 
are calculations, the dots are ex- 
perimental data from PHENIX 
collaboration [70] (the point for 
the central collisions is the upper 
limit for 90% C.L.). 
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Most of our results presented above are obtained considering a one unit 
rapidity window at midrapidity, while results for the full volume were pre- 
sented only for the SPS. Unlike the kinetic model of Thews et al. [73], our 
model does not contain dynamical aspects of the coalescence process. How- 
ever, in our approach, the width of the rapidity window does influence the 
results in the canonical regime. For the grand-canonical case, attained only 
at LHC energy, there is no dependence on the width of the rapidity window, 
due to a simple cancellation between the variation of the volume, proportional 
to the rapidity slice in case of a flat rapidity distribution, and the variation 
of iVcc"^! ^Iso proportional to the width of the rapidity slice. In Fig. 18 we 
present the centrality dependence of J/V' rapidity densities for RHIC energy 
and for different rapidity windows Ay from 0.5 to 3. The dependence on Ay 
resembles that of the kinetic model [73], but is less pronounced. The available 
data are not yet precise enough to rule out any of the scenarios considered. 
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However, for the kinetic model, the cases of sniall Ay seem to be ruled out 
by the present PHENIX data. We stress in this context that the size of the 
Ay window has a potentially large impact on the results at SPS energy. It is 
conceivable that no charm enhancement is needed to explain the data if one 
considers a sufficiently narrow rapidity window for the statistical hadroniza- 
tion. 



Particle Statistical pQCD NLO 
hadronization 



IF 0.228 0.155 

0.226 0.146 

D+ 0.190 0.095 

D7 0.189 0.089 

Ac 0.074 0.086 

Ac 0.074 0.062 



Table 4. Ratios of midrapidity 
densities for open charm hadrons 
relative to (D^ + D"), calculated for 
central collisions at LHC. The re- 
sults of the statistical hadroniza- 
tion model are compared to NLO 
pQCD calculations [72]. 




Fig. 19. Centrality dependence of rapidity densities of J/^p (per Npart) at LHC. 
Upper panel: sensitivity to dNch/dy, lower panel: sensitivity to dN^^^ /dy. 
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Ratios of mid-rapidity densities for open charm hadrons relative to (D*^ + 
D°) are presented in Table 4 for central collisions at LHC. The statistical 
hadronization model results are compared to NLO pQCD calculations [72]. 
In case of pQCD, the production of charm is identical to the elementary case, 
namely charm quark production in hard processes. Sizeable differences (up to 
a factor of 2, in case of Ds mesons) are seen. The measurements will certainly 
be able to distinguish between the two scenarios. 

In Fig. 19 we present the statistical hadronization model results on rapid- 
ity densities of J/tp per Npart for the LHC energy. We study the sensitivity 
on the two input parameters that are not well known at LHC, dNch/dy and 
dN^^^ /dy. Within the variations considered here (up to a factor 2 larger par- 
ticle multiplicities and a ±50% in the charm cross section) the changes in 
the yields are considerable, but the dependences on Npart remain the same, 
making this a rather solid prediction for LHC. 
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Fig. 20. Excitation func- 
tion of J/tp production in 
central nucleus-nucleus colli- 
sions. The symbols are sta- 
tistical hadronization model 
calculations. The lines are 
pQCD calculations for pp 
collisions [60], scaled for 
Npart=350, for two PDFs. 
The arrow denotes the exper- 
imental value [70]. 
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The excitation function of J/tp production (rapidity densities) is shown in 
Fig. 20. The statistical hadronization model results are compared to pQCD 
calculations for pp collisions [60], scaled for Npart=35Q, for two PDFs. Note 
that the PDFs, as well as the other inputs of the pQCD calculations [60] 
are different from those used to extract the charm cross section [62] which is 
an input to the statistical hadronization model. In any case, the yields are 
comparable in the two cases, implying that fine tuning of the input values 
will be needed to be able to distinguish between the two scenarios. As in the 
case of the total charm cross section [62] , the dependence of J/ ip production 
on the PDF choice is evident for the higher energies (LHC). 
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The results presented above were obtained under the assumption of sta- 
tistical hadronization of quarks and gluons. Wc have assumed that charm 
quarks are entirely produced via primary hard scattering and thermalized in 
the QGP. No secondary production of charm in the initial and final state 
was included in our calculations. Final state effects like nuclear absorption of 
J/V' [3] are also neglected. First RHIC data on J/ijj production support the 
current predictions, although the experimental errors are for the moment too 
large to allow firm conclusions. Also the RHIC results on open charm lend 
a strong support for this model. The statistical coalescence implies travel of 
charm quarks over significant distances e.g. in a QGP. If the model predictions 
will describe consistently precision data this would be a clear signal for the 
presence of a deconfined phase. We emphasize that the predictive power of 
this (and any similar) model relies heavily on the accuracy of the charm cross 
section, which is yet to be directly measured in nucleus-nucleus collisions. 

6 Outlook 

The field of ultrarelativistic nucleus-nucleus collisions has reached the stage of 
precision measurements, which arc able to provide fundamental information 
on the strongly interacting matter at high temperature and (energy) density 
and in particular on the quark-gluon plasma. It is now clear that such complex 
knowledge can only be achieved by a set of multi-faceted and complementary 
studies and that no one single observable is sufficient to characterize fully the 
properties of the QGP phase. Prom what we have briefly reviewed here it is 
clear that the global characterization of the collision has been convincingly 
achieved and, subject to further refinements, establishes beyond doubt that 
the conditions for the creation of QGP have been attained. We have shown 
that the study of particle ratios provide unique insight on the QCD phase 
diagram, while the study of charm hadrons provides a valuable glimpse into 
the QGP. 

We look forward to new high statistics and precision experiments from the 
SPS and, in particular, RHIC. From 2007 on the high temperature region of 
the phase diagram will be investigated in detail at the LHC with the dedicated 
ALICE experiment as well as within the ATLAS and CMS collaborations. 
The next serious attack on the high density-moderate temperature regime 
will be addressed further in the future with the new GSI accelerator facility. 
Interesting times are ahead! 
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